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^ ■ Abstract 

We have examined the properties of neutron-rich matter and finite nuclei in the 
^ I modified relativistic Hartree approximation for several values of the renormal- 



ization scale, fi, around the standard choice of fj, equal to the nucleon mass M. 
Observed neutron star masses do not effectively constrain the value of fi. However 
for finite nuclei the value fi/M=0.79, suggested by nuclear matter data, provides 
a good account of the bulk properties with a sigma mass of about 600 MeV. This 
value of fi/M renders the effective three and four body scalar self-couplings to be 
zero at 60% of equilibrium nuclear matter density, rather than in the vacuum. We 
have also found that the matter part of the exchange diagram has little impact on 
the bulk properties of neutron stars. 



1 Introduction 



The archetypal Walecka model is a point-particle, relativistic field theory in 
which nucleons interact via the exchange of a- and cj-mesons. If systems with 
non-zero isospin are studied the exchange of p-mesons also gives substantial 
contributions. This type of model has been widely employed to study matter 
both at normal and at high density e.g. [|1|. The Lagrangian takes the form 

C = N [i^d, - g^^u^ - ■ r + ^^'"duK ■r-M + g^a^ N 

+ \d,ad>'a - \mla'' , (1) 

where F^^ = d^u^ - d^u^, B^^ = df,K - d^h^ and af"" = ^ilYn"]- We 
have allowed for both vector and tensor coupling of the p field, b^, to the 
nucleon (the latter is only employed in evaluating the exchange or Fock di- 
agram). We note that the presence of the p meson renders the Lagrangian 
non-renormalizable, although a renormalizable model can be constructed [||]. 

In the mean field approximation, where the fields are replaced by their 
classical expectation value, the energy density for an infinite system is 

+ — [ dkk'^E* + — [ dkk'^E* . (2) 





Here the nucleon effective mass, M* = M — grjCr, E* = \/W+M*^ and 
Cf = {giM / rriiY , where i denotes the meson species. The subscripts N and 
P on the Fermi momentum kp and density n distinguish the neutron and 
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proton values. The total density n = tin + np, thus the neutron excess 
parameter a = ^^-^ = Normal nuclear matter corresponds to a = 0, 

but here we shall be interested in neutron star matter for which a ~ 1. Also 
matter with a ~ | is of interest for stellar collapse leading to supernova 
explosions. 

In the relativistic Hartree approximation (RHA), where the one loop di- 
agrams are taken into account, the energy density is the sum of the mean 
field contribution and a one-loop vacuum correction term £^rha = ^mfa + '^^- 
The one-loop vacuum correction term, caused by a shift in the single particle 
energies, receives contributions from the a-interaction only. It contains di- 
vergences and the necessary renormalization introduces a scale, fi. Explicitly 



A£{M*) = -iM*Un + Y, /„M^-"(M - M*)'^ . (3) 

V n=0 / 

d'-"Ag(M*) 



dM" 



In order to specify the finite parameters /„, we define \m{M*) - 
and require Ao(M) and Ai(M) to be zero so that the vacuum energy is zero 
and corresponds to a minimum. We can also set A2(M*) = at M* = M 
because different choices simply correspond to modifying the parameter 
which is fitted to equilibrium nuclear matter. The remaining parameters are 
specified by requiring that X^^M*) and A4(M*), the effective three and four 
body couplings, are zero at some scale and for simplicity we choose the same 
scale M* = in both cases. This yields 

/o = ln(/i/M), /i = l-41n(/i/M), 
/2 = -| + 61n(/x/M), /3 = f-4/i/M, /4 = -f • (4) 
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The vacuum correction term /S.£ can then be written in the form 
M / /i /i 



= f 1 - ^ + In ^ Vm - M*f 
7r2 V M MJ^ ' 



For the standard choice IH, /i = M, the first and second terms on the right 
vanish and one is left with the Chin-Walecka IHl result 



M*^ In — + M3(M - M*) - |M2(M - M 



*\2 



+ fM(M-M*)3-||(M-M*)^ . (6) 

Note that eq. (6) is designed to have no or contribution, such terms 
arise foifiy^M and are shown explicitly in eq. (5). We refer to this approach 
as the modified relativistic Hartree approximation (MRHA). 

In principle the results of an exact treatment would be independant of the 
renormalization scale, however this will not be the case when a simple ap- 
proximation such as the MRHA is employed. Therefore the renormalization 
scale should be chosen on physical grounds and we have previously found 

Q that in nuclear matter the choice /i 7^ M gives a better account of 
the second derivative of the energy/particle at saturation (the compression 
modulus) and also the third derivative (the skewness or anharmonicity) . Our 
purpose here is to test the MRHA further by applying it to finite nuclei and 
neutron stars. Nuclei are expected to provide a more stringent test, in partic- 
ular they are much more sensitive to the value of the nucleon effective mass 
than infinite systems. The value of M* in the interior of a heavy nucleus 
should approximate that found for equilibrium nuclear matter, which varies 
significantly with the renormalization scale. We will examine the predic- 



tions of the MRHA for the spherically-symmetric, doubly-closed-shell nuclei 
i6q^ 40Qg^ g^j^j^ 2°^Pb. In the case of neutron stars we are dealing with dense 
asymmetric matter. Apart from the obvious question of whether the MRHA 
supports neutron stars of the observed mass, we ask whether such global 
properties provide a significant constraint on the nuclear equation of state 
around the saturation density. It is also of interest to ask whether the pro- 
ton fractions are such that the direct Urea process, i.e. neutrino emission 
from a single nucleon, is possible. If so, this will provide the dominant cool- 
ing mechanism [^j since the standard modified Urea mechanism involves two 
nucleons. 

Exchange terms are a fundamental feature of nuclear physics, but they 
are not included in the MRHA. An additional goal of this work is to study 
the contribution of the lowest-order, two-loop exchange diagram to nuclear 
and neutron-rich matter. The vacuum contributions are large in magnitude, 
however it has been found that they become small if form factors are included 
to reflect the composite nature of the particles 0. We therefore consider 
just the matter, or Fermi sea, contributions which are much less affected by 
form factors. This is a tractable calculation which should allow us to get a 
qualitative indication of the importance of exchange effects in dense matter. 

The organization of the paper is as follows. Section 2 deals with the 
MRHA; we start with a brief discussion of nuclear matter in subsec. 2.1, 
turn to asymmetric matter in subsec. 2.2, neutron stars in subsec. 2.3. and, 
finally, finite nuclei in subsec. 2.4. In sec. 3 we investigate the effect of the 
matter part of the lowest-order exchange diagram for infinite systems. Our 
conclusions are given in sec. 4. 
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2 The Modified RHA 



2.1 Nuclear Matter 

For a given value of /i/M tlie parameters and are fitted to the equi- 
librium nuclear matter properties. We take a binding energy/particle of 16 
MeV. In the interior region of heavy nuclei the measured charge density in- 
dicates an equilibrium matter density of uq = 0.16 ± 0.01 fm~^; we choose 
the central value of 0.16 fm"'^. The parameter is fitted to the coefficient 
of the symmetry energy, a^y^ra = infj^ U=0! which is taken to be 30 MeV. 
The resulting parameters are listed in table 1 for a range of values of /x/M. 
The omega coupling constants, C^, agree qualitatively with the experimental 
result of 103 ±36 and are within the cited errors for /i/M= 0.73, 1.00 and 
1.25. Table 1 also gives the calculated compression moduli, K = kl^-^ {^^' 
However, Pearson has made the point that K is not uniquely determined 
by the limited amount of accurate breathing mode data available (see also 
Shlomo and Youngblood 0). He determines a linear correlation between K 
and the coefficient of the Coulomb term, KconiZ"^ / ^ in the leptodermous 
expansion. This coefficient can be conveniently written 



2 



^ coul r 

5ro 



s- 



(7) 



where tq is the radius parameter and S is the skewness or anharmonicity 
parameter defined by S* = kp-^ (^)' then has a linear correlation 

between S/ K and K as shown by the shaded band in fig. 1 (the box indicates 
the point obtained by Sharma [|l^, however changes little over the shaded 
area). The points give results ||^ obtained in the MRHA for various values 

6 



of the renormalization scale, The MRHA resuhs he within the error band 
for values of the renormahzation scale ft! 0.79 and 1.25. The other values of 
ji/M in table 1 lie off the shaded band and are included in our discussion for 
comparison purposes. The value 0.73 has the attractive feature that it is the 

self-consistent point where /i = M*^^. 

2.2 Asymmetric Matter 

We show in fig. 2 the binding energy/particle for nuclear matter and for pure 
neutron matter as a function of density; the full curves are obtained in the 
present MRHA approximation. At a given density the difference in energy 
between the pairs of curves is roughly independent of varying by less than 
10%. As regards the stiffness of the equation of state, it is often assumed 
that this can be deduced qualitatively from the compression modulus. We 
see, however, that it is not always a reliable indicator. Thus at high density 
the equation of state for ji — 0.73M is noticably stiffer than for ji — 1.25M, 
yet the compression modulus is smaller. A quantitative measure of the high 
density behavior of the equation of state is provided by the adiabaticity 
index, T = (dlnP/dlnn). At densities n > 2.5no the value of F for both 
nuclear and neutron matter lies in the range 2-3, with little dependence on 
the renormalization scale apart from the case /i — 0.73M. Here the pressure 
increases more rapidly with density and F is significantly larger than in the 
other cases. 

As we add neutrons to nuclear matter, i.e., as a increases from 0, the equi- 
librium binding energy is reduced and the corresponding density is reduced. 
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We have verified that this can be accurately described by the expression ITT 



risat(tt) = nnm(l - Ca ) . (8) 

The constant C is hsted in Table 2. Further one can define the compression 
modulus at the saturation density of the asymmetric system, i.e., the isobaric 
incompressibility Kiso, and to second order in a this can be written 



K,Ja) = K 



1 + {A-B)a^ =K l-Aa^ . (9) 



Here the coefficient A gives the asymmetry dependence of Kiso at nuclear 
matter density, while B allows for the change of equilibrium density in the 
neutron-rich system. These parameters, along with the combination A are 
given in table 2. Again we remark that this quadratic form provides a very 
accurate approximation, excepting only values of a close to the point where 
the system becomes unbound. The values of C and A are smallest for fi/M=l 
and show a factor of 3 variation with renormalization scale. Other calcula- 



tions [0 show a rather wide variation in the predicted values of C. However 
A is close to 1.6 in the chiral cr-model and it usually lies in the range 1.3-2.2 
in non-relativistic calculations, although larger values are to be found. These 
are comparable to the values in table 2. 

2.3 Stellar Matter 

We have calculated the properties of neutron stars using the hydrostatic 
equilibrium equations of Tolman, and Oppenheimer and Volkov [0. For 
n > 0.08 fm the MRHA equation of state is used. For densities 0.001 < 



n < 0.08 fm , we employ the equation of state of Negele and Vautherin [|T3], 



while for very low densities {n < 0.001 fm we use the Baym-Pethick- 



Sutherland [|15| equation of state. The calculation is carried out for charge 



neutral matter in beta equilibrium. The condition for beta equilibrium is 

dE{n,x) 

He = fJ-N - fJ-p = — , (10) 

where /ij, i = N,P and e, are the chemical potentials of the fermions, x = 
Z/A and E{n, x) is the energy per baryon. Muons will be present when /Zg > 
= 105.7 MeV, which is generally the case for n > hq. In the presence of 
muons, the proton fractions are determined by imposing the charge neutrality 
condition Ue + n^ = np and the energy conservation condition fif. = fJ-^- The 
lepton contributions to the energy density and pressure are given by Fermi 
gas expressions (the electromagnetic interactions of the leptons give negligible 
contributions) . 

The measured mass [|r^ of 4U0900-40, which is (1.85 ±O.3)M0, may pro- 
vide a limit for the maximum mass of a neutron star. The most accurate 
determination |T^, for one component of the binary pulsar PSR 1913 + 16, is 
(1.44 ± 0.01)Mq. The maximum masses in table 3 are consistent with these 
limits, although since hyperons are expected to reduce the maximum mass 



1^, [T^ the case with fx = 1.5M can probably be excluded and the value for 
fi = 1.25M appears dangerously low. Note that Mmax is not a monotonic 
function of K. Thus observed neutron star masses cannot be used to effec- 
tively constrain the compression modulus alone. Such properties of nuclear 
matter are not uniquely determined by the gross features of neutron stars 

20|] , which constrain the softness of the neutron-rich equation of state. In 



table 3 we also give the radius of the neutron star and the ratio of the cen- 



tral density to equilibrium nuclear matter density, the latter for a maximum 
mass star and a 1.44Mq star. The qualitative trends follow the high density 
stiffness of the equation of state. 

We turn now to estimate the maximum Keplerian frequency of rotation 



from the formula 121 



. 1 3 

M\ 2 / p \ - 2 

max I I ^ i-n * 



,,.^0,77^3=j lOV, (U) 

where Mmax and -Rmax are the maximum mass and corresponding radius for 
the spherical non- rotating star. The empirical relation (11) reproduces the 
results of more exact calculations ||22|, It is interesting to observe that 
changes in the renormalization scale serve to increase the Keplerian rotation 
frequency and give values which are on the upper end of those found with 
relativistic equations of state . 



Recently, it has been pointed out that the direct Urea processes 

n — > p + e~ + z7g , p + e~ — i> n + z/e (12) 

may constitute the principal avenues through which the late-time rapid cool- 
ing of neutron stars occurs. For these processes to occur in matter in which 
the only baryons are nucleons, momentum conservation requires that the 
magnitude of the electron concentration Xe exceeds a value |^ 

\xe\'/^>\{l-xy/'' -x'/^l, (13) 

where x = Z/A is the proton fraction. When leptons (e~ and fi~) are the 
only source of negatively charged particles, charge neutrality and eq. (|1^) 
stipulate that x exceed a value in the range 0.11 — 0.15. In the models 
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considered above, the symmetry energy rises sufficiently rapidly with density 
so that the calculated proton fractions are large enough for the direct Urea 
cooling mechanism to be operative. 

Since the calculated maximum neutron star masses exceed the observed 
values for all but very large values of /i/M, we conclude that the gross prop- 
erties of neutron stars do not effectively constrain the renormalization scale 
H/M. 



2.4 Finite Nuclei 
2.4.1 Formalism 

For finite nuclei the Coulomb interaction is needed so that the Lagrangian 
(1) is supplemented by the photon Lagrangian 

Cph = -lUuF" - eiVi(l + Tsh'-A^N . (14) 

Here is the Maxwell field and /^i, is the electromagnetic field strength 
tensor. It is straightforward ^ to obtain from C + Cph the equations of 
motion in the mean field approximation. The vacuum is taken into account 
by means of the local density approximation and the leading correction terms 
of the derivative expansion [E^ . The latter give small, although not negligible 



corrections p6[. The derivative expansion appears to be rapidly convergent 
[p^ , [28[] so that the lowest order corrections should suffice. 

The corrections of the derivative expansion can be obtained by evaluating 
the one-loop polarization contribution to the scalar and vector meson prop- 
agators [^. The question arises as to where these propagators should be 
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renormalized, at M* = M or at M* = /i. If in the former case one chooses 
A2(M) = and in the latter A2(/u) = 0, where A2 is defined in sec. 1, then 
changing the renormahzation point for the scalar meson is simply equiva- 
lent to varying m^, which is a free parameter anyway. For the u;- meson it 
is reasonable to require that, in vacuum, the pole of the propagator lies at 
the known physical mass. Then changes in the renormahzation point will 
not affect the predictions, but rather the interpretation of the parameters. 
Therefore it is sufficient to take the derivative contributions in standard form 



[ p5|l so that the complete scalar density is 
(2ja + 1: 



GJr) 



FJr) 



+ 



dA£ 



dM* 



An' 



2g In — (— + - — 

M \ dr"^ r dr 



9^ 
M 



2 ^da\_ 2gl (du^ 



dr I 3M* \dr 



(15) 



where G and F are the components of the Dirac spinors for the occupied 
states and M* is the spatially dependant effective mass, M*(r) = M —g„a{r). 
Only the timelike component of uj^ is non-zero for the spherically symmetric 
nuclei we consider. The full vector density is 

(2j„ + 1) 



37r2 



GairY + F^ir) 



M* ( d'^ujQ 2 dujQ \ guQu) ( da dujQ 



dr"^ ^ r dr 



M* 



dr dr 



(16) 



and the contribution to the total energy from the vacuum is given by 



AE = in r^dr 



A£{M*) 



A. 



M* (daX gl , M* ( dujQ 



In — + ^In 



M \ dr 



o In- 

67r2 M 



dr 



where M* = M*{r). 
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2.4.2 Results 

We now consider the predictions of the MRHA for finite nucleiH, specifically 
the doubly closed-shell nuclei ^^O, '^''Ca and ^°^Pb. In addition to the param- 
eters of table 1, the meson masses are required. We use the physical masses 
for the fh and tu-mesons; the value of mo-, which is not a priori known, is 
discussed below. A pure vector pN coupling is employed in these calcula- 
tions. As we have remarked, the vacuum is treated by means of the local 
density approximation supplemented by the leading terms of the derivative 
expansion. The effect of the latter are fairly small, giving a decrease in the 
binding energy, an increase in the radius for a fixed a mass and very little 
change in the single particle energies. If mo- is adjusted so that similar radii 
are obtained with and without the derivative terms, the net result is a small 



increase in the binding energy when the derivative terms are included \2t 



We shall present results which include the contributions of the derivative 
terms. 

In table 4 we summarize the bulk properties of the nuclei as a function 
of fi/M for a fixed value of mo= 550 MeV. The theoretical values of the 
binding energy/particle include a correction for the cm. kinetic energy 
The charge densities, and therefore the radii quoted, are corrected for the 



finite size of the proton and for cm. effects ||3^. In table 4 we also give 
the separation between the neutron 2/ levels in ^°'^Pb so as to indicate the 
order of magnitude of the spin-orbit sphttings. In studying table 4 it must 



"'^This is also under investigation by Blundcn |30| | with similar results where comparison 
can be made. 
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be borne in mind that while is fixed from nuclear matter, the value of 
m„ can be regarded as a free parameter. Decreasing m„ has the effect of 



increasing and decreasing and vice versa p4|. Thus for jj/M =0.73 



the charge radius of Pb is too large, while those of O and Ca are too small 
so that we shall not be able to obtain a reasonable account of these nuclei 
by varying rrif^. For yu/M=1.5, the deviations of the radius from experiment 
are quite disparate for Ca and Pb and this cannot be resolved by adjusting 
m„. Further in this case the spin-orbit splittings are much too small- the 
value in table 4 is less than a third of the experimental figure. This is due to 
the large value of the effective mass at saturation in nuclear matter (0.86). 
As is well known [|l| smaller effective masses ~ 0.6 yield better values for 
the spin-orbit splittings and this is evident from comparison of tables 1 and 
4. The value /i/M=1.25 also yields quite small spin-orbit splittings due to 
a large value of M*^^, although, for comparison purposes, we will consider it 
along with the values 0.79 and 1.0 in subsequent calculations. 

The results obtained by varying so as to approximately fit the radius 
of Ca and, if possible, also Pb are given in table 5. The binding energies for 
^/M=0.79 are in strikingly good agreement with the data for a reasonable 
value of mo-=600 MeV. With yu/M=1.25 the binding energies are low and 
the required a mass is unreasonably small. The value /i/M=1.0 also yields 
rather low binding energies and further one cannot obtain a good account of 
the charge radii of Ca and Pb simultaneously. We should emphasize at this 
juncture that our results are obtained by fitting to a nuclear matter saturation 
density, uq = 0.16 fm~^. A number of other authors, e.g. [|l], ^ have 



used the value 0.15 fm ^ which represents a reduction in kp of 2%. This 
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has the effect of increasing the nuclear radius by approximately 2% so that a 
larger mg. can be used and this increases the binding energies. For example, 



Wasson takes 550 MeV {j2/M=1.0) which gives binding energies 1-2 
MeV short of the experimental value and in fact for /i/M=0.79 the nuclei 
are overbound with this saturation density. In both cases close agreement 
with the observed radii for all three nuclei cannot be obtained. 

Adopting the values of m^j in table 5 which are appropriate to uq = 0.16 
fm~^, we obtain the single particle energies of ^°^Pb which are shown in figs. 
3 and 4 for neutrons and protons respectively. We indicate both occupied and 
unoccupied experimental states and their theoretical counterparts, if bound, 
for levels near the Fermi energy. The shell closure and rough level ordering 
are reproduced quite well, however the occupied neutron (proton) levels are 
too strongly (weakly) bound which could be improved by increasing the p 



coupling p4|. For example an increase in of 40%, corresponding to a 
symmetry energy of 35 MeV, would shift the occupied single particle levels 
in the right direction by about 1 MeV; the unoccupied levels are shifted by 
a somewhat smaller amount. 

The spin orbit splittings that we find are too small; on average they 
are 80, 60 and 30% of the experimental values for /i/M=0.79, 1.0 and 1.25, 
respectively, and the percentage becomes somewhat smaller for the lighter 



nuclei. For the standard renormalization. Fox ||3^ has made the same point 
and noted that the mean field results are close to, or a little larger than, 
the data (see also, for example, the fits in ref. [^). Finally we show in 
figs. 5-7 the comparison of the predicted charge distributions with the data 
(experimental errors are not indicated since they are negligible except near 
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the center of the nuclei). In Pb for yu/M=1.0 our central density is too high, 
whereas Fox |Q achieves better agreement, presumably due to his choice of 
hq = 0.15 fm~^. We do however find a region of positive slope, although 
it not as pronounced as the data; mean field calculations do better in this 
respect. Our results for /i/M=0.79 and 1.25 are significantly better, both for 
radial distances of 2-5 fm and in the important tail region. For Ca, /i/M=1.0 
clearly gives the best agreement with the data, although the other cases are 
not unreasonable. In O the different theoretical results are of comparable 
quality. 

Viewing the results for finite nuclei overall, the renormalization scale 
fi/M=0.79 gives the best representation of the data. The additional terms 
which are introduced into the vacuum correction AS were given in eq. (5). 
This value of /i/M renders the self-consistently determined effective three and 
four body scalar self-couplings to be zero at a nuclear density of n = 0.092 
fm^^, which is roughly halfway between the equilibrium point and the vac- 
uum. With fi/M =0.79 the coefficient of the induced cr^ term is negative and 
this has long been known to be necessary in mean field treatments which fit 



parameters of a and a terms to nuclei [34, 35, 36 
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3 Exchange Contributions for Infinite 
Systems 



3.1 Formalism 

We turn now to assess the importance of exchange contributions from the 
two- loop diagram shown in fig. 8; here a sohd hne represents a nucleon with 
propagator S and a dashed hne represents a meson with propagator D. Along 
with the a and u mesons we shall include the p since this is clearly necessary 
for neutron-rich matter. The exchange diagram gives a contribution to the 
energy density 

where k = pi — p2 and Lorentz indices have been suppressed. We use the 
mean field form of the propagators, treating the two-loop contribution as 
a perturbation. The matter part of this contribution to the energy density 
involves the basic integrals 

1 

J(w;,u,6) = i dv(l--)\n^ ^- . 19 

J \ J [u ~ vy + uvw 

For completeness, we first give the expressions for a- and cj-exchange. 
The energy density is the sum of neutron and proton contributions 

£u; = (kpN) + (kFp) ■ (20) 
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Here 



9l 



647r4 

,2 



327r4 



where 



E 



*2 



= (m^/M*)^ and = (m^/M*)^ 



The energy density for p exchange can be written 



Sp = Cp {kpN, kpN) + Cp {kpp, kpF) + 4ep (/cfat, /cfp) . 



(22) 



For the general case with both vector and tensor couphng each Cp is the sum 
of vector- vector, tensor-tensor and vector-tensor contributions 

ep {kpi, kpT) = e'p {kpi, kpT) + e" {kpi, kpT) + e^* {kpi, kF2) ■ (23) 
These contributions can be written in the form 

.2 1 



e7 {kFi,kF2) 



9p 

2 y 327r4 



(A;fi£;>i - M*2 In ^i) (A;p2£^F2 - M*^ In 6) 
-M*^(2 + ^,)7K,6,6)1 , 



fp \ 



^4My 647r4 

r*2 



;iO + Wp)(A;nE^i-M*Mnei 



x{kp2EU - M*' In 6) - - M*'wp{8 + Wp)I{wp, 6, 6) 



e;'{kFi,kF2) 



8M / 167r4 



(A;pi£;;i - M*' \nC,){kF2E*F2 - M*' In 6) 



M*X/K,ei,e2) 



(24) 
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where Wp = {nip/M*)'^ and the subscripts 1 and 2 on Ep and ^ refer to the ap- 
propriate Fermi momenta. If the no-retardation approximation is made, i.e., 
energy differences {EI — E2) are set to zero, this can be put in the form given 
by Bouyssy et al. [|12|. Note that they remove the (three-dimensional) delta 
function contribution from the tensor-tensor piece, arguing that it should be 
suppressed by short-range correlations. 

Thus in this approximation the total energy density £^totai = ^^rha + ^^exch, 
where £^exch = £u + £lu + £p- For a given density the value of M* is obtained 
from the usual minimization condition dStotai/dM* = and here both the 
RHA and exchange contributions are taken into account. A word of caution 
is in order regarding the high density behavior, kp 00. Since in this 
limit £^MFA oc ra^, while £a and oc ||^ these exchange contributions 
are negligible. However the p tensor coupling involves a derivative which 
brings in an extra power of the momentum. This causes the tensor-tensor 
contribution to dominate the asymptotic behavior so that Sp oc —n^ and the 
total energy density £ oc (C^ — ^^~g ^ fp)n'^ has the possibility of becoming 
negative for sufficiently large fp. Furthermore the equation for the effective 
mass M* will ultimately be dominated by the p vector-tensor and tensor- 
tensor contributions which will give rise to a negative value of M* (for positive 
fpl dp)- Such problems arise in principle for the p, but only for densities > 
15^0 which are well beyond our range of interest and probably beyond the 
range of applicability of a simple baryon-meson picture. 
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The Landau effective mass, M£, is defined in terms of the velocity at the 
Fermi surface, vp, and the single particle spectrum, e(A;), according to |^ 



Vp 



where vp 



de{k) 



dk 



k=kf 



(25) 



The well known result Ml = {kp + M*'^)^ is obtained in the MRHA. Here we 
shall study the modification of induced by the exchange diagram for the 
case of nuclear matter where the neutron and proton Fermi momenta are the 
same. In order to perform the calculation we need e{k). This can conveniently 
be obtained by taking the functional derivative of the energy density with 
respect to the occupation probability, 6S/6n(k), where, for nuclear matter, 
n(k) = {A7r^)^^6{kp — |k|). The contributions of scalar and vector exchange 
to e{k) have been given by Chin [Q]. The contribution of the p can be written 
as the sum of vector-vector, tensor-tensor and vector-tensor terms. 



e^{k) = ey{k)+el\k)+e;\k). 

For nuclear matter these terms are 
>2 3 



(26) 



C(fc) 



9p _ 

2 J Stt^E; 



M 



*3 



AM) Stt^E^ 



kpE*p - M*' In e - (2 + «;,) J {wp, ^fc, 



24^* - 3M*2(5 + lwp){kpE*p - M*2 InO 



+ — ; Wp{l + \wp)J {wp,^k,0 



8M J Att^E; 



kpE*p - M*^ In^ 



*3 



k 



-WpJ iWp,^k,0 



(27) 



where Ef = k^ + M*^, = {k + El)/M* and the other quantities are defined 
above. The derivative of e is then easily evaluated to obtain the contribution 
to the Fermi velocity. 
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3.2 Results 

When the exchange diagram is included it is necessary to refit nuclear matter 
properties. In order to restrict the number of parameters, we take a a mass of 
550 MeV and fix the p-nucleon couplings. For the vector part we take 5'p/47r = 
2.4 (Cp = 44.9) from vrA^ scattering [Q, assuming p universality. Reference 
indicates a tensor- vector ratio, fp/gp, of 6.6, which is at the upper end 
of currently accepted values 0. We shall also consider the smaller ratio of 
3.7 derived from the isovector anomalous magnetic moment of the nucleon in 
the vector dominance model; studies of nuclear systems in the Hartree-Fock 



approximation favor the smaller value. The a and uj coupling constants 
which result from fitting the saturation properties are collected in table 6. 
Since the net effect of exchange is repulsive, it is necessary to reduce the u) 
coupling and for /i/M=1.5, is unreasonably small and, in fact, we were 
not able to obtain a satisfactory fit with the stronger tensor coupling. We also 
show in table 6 the compression moduli which show some modest changes 
from their previous values; for /i/M = 0.79 and 1.0, K is reduced, whereas 
in the other cases K increases. In the last column we give the coefficient 
of the symmetry energy, asymm- Here the exchange contribution from the a 
{(jj) is positive (negative), while the effect of the p is quite small; the net 
contribution from exchange is large, ~ +15 MeV, as is well known. Thus 
without the direct p contribution flsymm ~ 30 MeV and, adding in the p, one 
obtains a value in the neighborhood of 40 MeV. Given that currently accepted 



values are in the range 27-35 MeV [^, our results are somewhat high 
particularly for p/M < 1. 
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The effect of exchange on the binding energy/particle as a function of 
density is shown in fig. 2. The dashed curves correspond to the weaker 
tensor p couphng for which the effect of exchange is fairly small, although the 
impact of the increased symmetry energy is noticeable for neutron matter. 
The dash-dotted curves show the results for the stronger tensor coupling; 
the reduced values of here lead to a softening of the equation of state 
for densities beyond 0.43 fm~^ in all cases. The parameters C and A of eqs. 
(8) and (9) characterise the change in equilibrium density and compression 
modulus with the neutron excess parameter a. The effect of exchange upon 
these quantities is case dependant, but typically the values of C in table 2 
are increased by 10% and the values of A are increased by 20%. 

In table 7 we show the individual contributions of the MRHA and ex- 
change to the Fermi velocity (see eq. (22)). The exchange contribution arises 
from (T (negative) and ui and p (positive) contributions, none of which is in- 
significant. The total exchange contribution is about 5% (25%) for the weak 
(strong) tensor coupling. In the weak case the net results for the Landau 
effective mass are quite close to those given in table 1 for the pure MRHA, 
whereas a reduction of ~10% is obtained with the strong tensor coupling. In 
neither case are the exchange contributions to the effective mass sufficiently 
small that they can safely be ignored. 

Turning to neutron stars, we find that properties given in table 3 are 
little changed when exchange effects are included using the weaker tensor p 
couphng. With the stronger tensor couphng the equation of state is softer 
which leads to a modest increase in the central densities and a small decrease 
in the maximum masses. The predictions for p/M=0.79 and 1.0 become very 
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similar, in particular the maximum mass is close to 2.2QMq. The predictions 
for yu/M=0.73 and 1.25 are also similar and the maximum mass of 1.79Mq is 
probably a little low in view of the expected effect of hyperons . The 

cooling of neutron stars by the direct Urea process is allowed in all cases, as 
was the case without exchange. 

Summarizing this section, we have examined the matter contribution aris- 
ing from the lowest order exchange diagram using the minimum number of 
mesons and found quite modest effects for neutron-rich systems, even with 
the strong tensor- vector ratio for the p coupling. It is unlikely that inclusion 
of, for example, the pion field would qualitatively change this conclusion. In- 
deed rather small effects were obtained from the Fock terms in (pure matter) 
Hartree-Fock calculations which included the pion as well as a number 
of baryon fields. 



4 Conclusions 

The main thrust of our work has been to examine the implications of different 
choices of the renormalization scale in neutron-rich matter and finite nuclei, 
following our previous work [|2|, Q with nuclear matter. Physical results would 
be independant of this scale if we had embarked on a complete renormaliza- 
tion program which accounted for all the important physical effects, but, of 
course, we cannot claim that this is the case. Within the confines of the 
relativistic Hartree approximation we have previously found that nuclear 
matter favors values of fi/M=0.79 and 1.25. The observed masses of neutron 
stars do not distinguish different values of p/M; however, for finite nuclei 
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we find that the phenomenology with yu/M=0.79 is favored since it provides 
rather reasonable agreement with the wide range of data considered. We 
can compare with the mean field approach (ref. and references therein) 
which has the disadvantage that it excludes the vacuum which is an inherent 
feature of relativistic field theories. Provided and cr^ terms are included 
in the Lagrangian, this approach gives somewhat better agreement for finite 
nuclei, particularly as regards the single particle energies. However it must 
be borne in mind that non-relativistic calculations |^4| show that correlations 
play an important role in determining the spectral strength distributions in 
actual nuclei. 

Changing the renormalization scale /i/M from 1.0 to 0.79 can be viewed 
as choosing the effective three and four body scalar self-couplings to be zero 
at a point which correspond to 60% of equilibrium nuclear matter density, 
rather than in the vacuum. For /x/M< 1, the coefficient of the a'^ term is 
negative and one can worry that the effective potential of the scalar field, 

[/eff(a) = Im'y + AS{M- g„a) , 

has quite a shallow minimum 0] at the vacuum value of cr = for yu/M=0.79. 
More troublesome is the fact that for a ^ — oo the effective potential is un- 
bounded from below, f/gfj — oo, irrespective of the choice of the renormal- 
ization scale. This means that formally we are dealing with a local minimum 
rather than an absolute minimum in the energy. Thus, while the phenomenol- 
ogy is successful in explaining a large amount of data, unresolved theoretical 
issues remain. 

We thank C.J. Horowitz for giving us a copy of his computer code for 
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computing time from the Minnesota Supercomputer Institute is gratefully 
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Table 1 

MRHA coupling constants and saturation properties 



M 


M 


K 

(MeV) 


Cl 


Cl 


Cl 


0.73 


0.73 


157 


132.8 


279.4 


82.5 


0.79 


0.66 


354 


180.6 


317.5 


73.5 


1.00 


0.73 


461 


137.7 


215.0 


81.6 


1.25 


0.82 


264 


78.6 


178.6 


90.8 


1.50 


0.86 


189 


51.2 


175.6 


94.4 



Table 2 

Parameters specifying the saturation density and isobaric 
compression modulus of neutron-rich matter 





C 


A 


B 


A 


M 








0.73 


1.62 


3.36 


5.45 


2.09 


0.79 


0.76 


1.56 


4.45 


2.89 


1.00 


0.54 


1.05 


2.16 


1.10 


1.25 


0.91 


1.75 


2.96 


1.21 


1.50 


1.25 


2.39 


3.64 


1.25 
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Table 3 

Neutron star properties in the MRHA 





K 




R 


ric 


ric 




M 


Me 


no 


no 




(MeV) 




(km) 




(I.44M0) 


(10^ s-i) 


0.73 


157 


2.18 


11.2 


6.2 


2.6 


0.96 


0.79 


354 


2.53 


12.7 


4.8 


1.9 


0.86 


1.00 


461 


2.30 


12.1 


5.3 


2.0 


0.88 


1.25 


264 


1.86 


10.7 


7.3 


3.1 


0.95 


1.50 


189 


1.59 


9.9 


9.1 


4.5 


0.99 



Table 4 

Bulk properties of nuclei for rua — 550 MeV 



M 




BE 
A 

(MeV) 


Teh 

fm 


Ca 

BE 
A 

(MeV) 


Teh 

fm 


BE 
A 

(MeV) 


Pb 

Teh 

fm 


2f u split 
(MeV) 


0.73 


9.42 


2.66 


9.52 


3.42 


8.33 


5.57 


1.09 


0.79 


7.52 


2.80 


7.99 


3.54 


7.47 


5.54 


1.64 


1.00 


6.08 


2.64 


6.92 


3.38 


6.70 


5.40 


1.31 


1.25 


9.16 


2.51 


9.26 


3.28 


7.99 


5.40 


0.73 


1.50 


10.59 


2. 16 


10.31 


3.21 


8.56 


5.13 


0.50 


Expt. 


7.98 


2.73 


8.45 


3.48 


7.86 


5.50 


1.77 
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Table 5 



Bulk properties of nuclei with fitted 



(MeV) 




— r-ch 

(MeV) fm 


Ca 

— r-ch 

(MeV) fm 


Pb 

— r-ch 

(MeV) fm 


0.79 600 
1.00 450 
1.25 350 


8.13 2.74 
4.03 2.79 
5.53 2.77 


8.47 3.49 
5.35 3.49 
6.65 3.48 


7.73 5.51 
5.82 5.45 
6.63 5.50 


Experiment 


7.98 2.73 


8.45 3.48 


7.86 5.50 



Table 6 

MRHA+ exchange coupling constants and saturation properties 



M 


K 
(MeV) 


CI 


CI 


9p 


'j^symm 

(MeV) 


0.73 


191 


88.3 


257.1 


3.7 


41.7 




204 


57.9 


212.0 


6.6 


38.6 


0.79 


332 


128.9 


283.6 


3.7 


42.2 




274 


96.4 


247.1 


6.6 


39.7 


1.00 


460 


106.7 


201.5 


3.7 


34.8 




421 


80.9 


186.6 


6.6 


34.4 


1.25 


275 


49.5 


169.6 


3.7 


35.4 




280 


26.3 


149.5 


6.6 


34.6 


1.50 


198 


18.7 


169.3 


3.7 


37.9 
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Table 7 



Landau effective masses in the MRHA+exchange approximation 







fplOp = 3.7 






fplgp = 6.6 






Contribution to vp 


Ml/M 


Contribution to vp 


Ml/M 




MRHA 


exchange 




MRHA 


exchange 




0.73 


0.350 


0.011 


0.78 


0.339 


0.086 


0.66 


0.79 


0.378 


0.021 


0.70 


0.364 


0.089 


0.62 


1.00 


0.354 


0.038 


0.72 


0.348 


0.103 


0.62 


1.25 


0.322 


0.023 


0.81 


0.318 


0.095 


0.68 


1.50 


0.308 


0.012 


0.87 
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Figure Captions 



Fig. 1. The ratio of tlie tliird and second derivatives of the nuclear matter 
equation of state at the equilibrium point, S/K, versus the compression 
modulus, K. The shaded band is obtained from the breathing mode data 
[§] and the filled square is the value suggested by Sharma |]lOl- The MRHA 
points are labelled by the value of the renormalization scale, [ijM. 
Fig. 2. The binding energy/particle as a function of density for neutron and 
nuclear matter with four values of the renormalization scale. The full curves 
correspond to the MRHA and the dashed (dash-dotted) curves correspond 
to MRHA + exchange with the weaker (stronger) tensor p coupling. The 
upper (lower) set of curves are for neutron (nuclear) matter. 
Fig. 3. Occupied and unoccupied neutron levels near the Fermi energy in 
208p|^ The experimental data are compared with predictions for yu/M=0.79, 
1.0, 1.25. 

Fig. 4. As for fig. 5, but for protons. 

Fig. 5. Comparison of the experimental charge density for ^^O with 

theoretical predictions for /i/M=0.79, 1.0, 1.25. 

Fig. 6 As for fig. 5, but for ^°Ca. The data are from ref. p8 |. 

Fig. 7 As for fig. 6, but for ^ospb. 

Fig. 8. The exchange or Fock diagram. In our calculations the dashed line 
represents the exchange of a cr-, a;- or p-meson. 



33 



